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FOREWORD 


1. REVIEW AND PROBLEM INDEX 


Project Whirlwind makes use of the facilities of the Digi¬ 
tal Computer Laboratory. The principal objective of the Project is 
the application of an electronic digital computer of lar-e capacity 
and very high speed (Whirlwind I) to problems In mathematics, science 
engineering, simulation and control. 

Whirlwind I Is of the high-speed electronic digital type, 
in which quantities are represented as discrete numbers, and complex 
prohloms are solved by the repeated use of fundamental arithmetic 
and logical (l.e., control or selection) operations. Computations 
are executed by fractional-microsecond pulses In electronic circuits 
of which the principal ones are (!) the flip-flop, a circuit cont.ln- 
Ing two vacuum tubes so connected that one tube or the other Is con 
luctlng, but not both, ( 2 ) the gate or coincidence circuit; ( 3 ) the 

:r ro t 0, ln which ,,in “ ry dieua ’ ,or,-d »» 

two directions of magnetic flux within ferro-magnetlc cores. 

, a . e " n ' ,rl, " nd 1 U9BS numbers of 16 binary digits (equivalent 
to about 3 decimal digits). This length was selec'ed to Jl.lt th^ 
m ch n„ to a practical slxo. bu, It permits the computation of Tny 
simulation problems. Calculations requiring greater number length 
are handled by the use of multiple-length numbers Kapld-access 
magnetic-core memory has a capacity of 32,768 binary digits Pre- 

ToT on a r a v : hc coBputor 18 • ,o ' , ’°° "‘"«>«-‘- dd -„ v ope g ri Ls p: r 

second, equivalent to about 20.000 mu. t Ipl let ions per second. 


23. 1936. ^InTm, 0 ^™ °V1936 to December 

hours of Whlrwlnd computer time allocated to 399,46 hour * of the W9W.06 

tat Iona (8*EC) Croup Of the S94.06 al^c«J “ d En « ln - r ‘"« Compu- 

tlm. becauae of computer malfunction. 9 J^ho^a co - pu ‘T . 2.81 ... down 

for terminal equipment te.ting and calibration dll »lloc.t«l time were u.ed 

run operations not logged to Zpecmc pJo" ^a - ” t ' ,tl0 “ “ d *•*■“»• ‘"ter 

Th. r..ui,r , of 6 2o p orj^r*p^iiL ri ^vr^ u or"n o ^ r i'T 18 ° r 

In these 20 problems, there are repre.ented 17 b ®, lnclu<led ln •cad.mlc these. 

Bachelor .. Twenty-two of the ornhi... i - ” Uoctor * 1 theses, 5 Master a and 1 
•ored at MIT by th. Office of Naval R...."h. 0rlgl “ te<l fr °“ rMMrch ***•«*• «Pon 

ports have been'si^n te/' ‘ J - ”.' 0 th"* Pr ° ble " ,or whlch Progr... re- 

to th. f,., d of appllcmtloo and ,h.“"* !'1. Pr °. b1 '" ‘ r ° * rr,,n « < * d ‘wording 

«n Tab.. 2-„. th . problems “ ll.t^acco^n^tn t ‘"* U88d "" «!*«■ 

blem involved ln each l n each table th. >.!?* 1 , h * P rl »c‘P«l mathematical pro- 
wheth.r the problem 1 . for ac.domlc credit and "h Tf'" tbe probl “ "“"Bor indicate, 
explained ln .ectlon 21 , Introduction" -1 * ^ ’"“h" “ “ ■ p “' d " d The cod. ,. 
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WHIRLWIND CODING AND APPLICATIONS 


2. WHIRLWIND CODING AND APPLICATIONS 


1.1 Introduction 

TUB letters have the following significance. 


106 C. NIT SK1SN1C PROJECT 

The relatively little coaputatlonal work we have done this last quarter con¬ 
sisted of cross correlation on a selsalc record to locate an echo position; testing 
the stationarlty hypothesis on several selsalc traces and selsalc well logs; and be¬ 
ginning the construction of theoretical randoa elastic aedla with desired scattering 
properties. Work will be continued and expanded In the next quarter. 

Prograamers Working on this problea are S. Treltel and D. Grins. 

S. Stepson 

Geology and Geophysics 


120 B,N. THERMODYNAMIC AND DYNAMIC EFFECTS OP WATER INJECTION INTO HI0H- 
TEMPERATURB, HIGH-VELOCITY GAS STREAMS 


A laplles the problea Is NOT for acadealc credit, Is UNsponsored. 

1 ) laplles the problea IS for acadealc credit. Is UNsponsored. 

C lapllos the problea is NOT for acadealc credit, IS sponsored. 

D laplles the problea IS for acadealc credit, IS sponsored. 

N laplles the problea la sponsored by the Office of Naval Reaearch. 

L laplles the problea Is sponsored by Lincoln Laboratory, 

The absence of a letter Indicates that tho problea originated within the SSJSC 
Group. 


2.2 ProgroiH Reports 


100 COMPREHENSIVE SYSTEM OP SERVICE ROUTINES 

During the past quarter, the Coaprehenslve Systea has been running at Its 
usual routine pace. 


Digital Coaputer Laboratory Staff 


This problea Is connected with the developaent of a potential gas turbine 
component, called an aerotheraopressor', in which a net rise In stagnation pressure 
of a hot gas stream Is brought about by cooling of the stream. The cooling of the hot 
gas Is accomplished by evaporation of liquid water which is Injected Into a high- 
velocity region of the flow. The concepts underlying the operation of the aorotherao 
pressor are an outgrowth of comparatively recent work in the field of gaa dynaaics, 
and its proposed function In the gas turbine cycle Is analogous to that of the con¬ 
denser In a steam power plant. 

The device consists of a converging nozzle which accelerates the exhaust 
gases from the turbine into a circular duct of varying diameter terminated by a con¬ 
ventional conical diffuser, which recovers the kinetic energy of the flow before dls 
charging It to the atmosphere. At the entrance of the duct, special injectors deliver 
minute Jets of water which are In turn atomized by the rapidly moving gas stream. 

The changes In state within the aerotheraopressor are brought about by the 
simultaneous thermodynamic and dynamic effects of (a) evaporation of the liquid water, 
(b) momentum and energy Interactions between the phases, (c) friction, and (d) varia¬ 
tions In cross-sectional area of the duct. Under proper circumstances, these effects 
bring about a not rise In stagnation pressure across the device, further descriptions 
of this device may be found In earlier reports, beginning with Summary Report No. 32, 
Fourth garter, 1952. 

Whirlwind I has played a very Important role in the development of the aero¬ 
theraopressor. Tho theoretical one-dlmenslonal analysis of the aerotheraopressor 
process Involves the simultaneous solution of seven non-linear first-order differential 
equations for a wide variation of initial conditions, Coaputatlon of results In a 
reasonable time by any method other than a fast digital machine such as Whirlwind 
would have boon impractical. 
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The agreement obtain* between experimental and theoretical re.ult. 1. good. 
Detailed compartaona are given In the following report*: 

1) Gevrll B. D. , A Theoretical InventlgatIon of the Thermodynamic and 
Dynamic Effect, of Water Injection Into High-Velocity, High-Temperature Ga. Stream.. 


2) Powle, A. A., An Experimental Inve.tlgatlon of an Aerothermoproa.or 
Having a Ga. Plow of 25 Pound. Per Second." 

3 ) Erickson, A. J., "A Theoretical and Experimental Inve.tlgatlon of the 
Aerothermopre.aor Proceaa." 


The experimental and theoretical work on the Aerothermopre.aor Project wa. 
carried out at M.I.T. under the .pon.or.hlp of the Office of Naval Ke.earch, and ... 
directed by Profe..or A. H, Shapiro of the Department of Mechanical Engineering. The 
theoretical aepecta of the problem treated by Whirlwind 1 were programed by D. D. 
Oavril and carried out by Gevrll and Krickeon. 


A. J. Erlckaon 
Mechanical Engineering 


126 C. DATA REDUCTION 

Problem 126 la a very large data-reductlon program for uao In the Servo¬ 
mechanisms Laboratory. The over-.ll problem la compo.ed of manv component section, 
which havo been developed separately and are now being combined Into complete proto¬ 
type programs. Description, of the various component soctions have appeared In past 
quarterly report*. After the development and teatlng of the prototype Whirlwind pro¬ 
grams la completed, the programs will be recoded for other, commercially available, 
largo scale computers (probably the ERA 1103, IBM 701 and IBM 704 computers), for use 
by Interested agencies for actual data reduction at other location*. The program, are 
currently being developed by Douglas T. Ross. David F McAvlnn, Benson H. Scheff, and 
Dorothy A. Hamilton, Servomechanism* Laboratory staff members, with the assistance of 
John F. Walsh. This work Is sponsored by the Air Force Weapon* Guidance Laboratory 
through Division of Sponsored Research Project 7138, 

The nature of the problem requlros extreme automatlclty and efficiency In 
the actual running of the program, but also requires the preaence of human operators 
In the computation loop for the purpose of decision making and program modification. 
For this reason, extensive use Is made of output oscilloscopes so that the computer 
can communicate with the human, and manual intervention registers so that the human 
can communicate with the computer In terms of broad Ideas, while the computer Is 
running, and have the comouter program translate these Ideas Into the detailed steps 
necoasary for program modification to conform to tho human operator's decision. The 
program which does this translation and aodlflcatlon Is called the Manual Intervention 
Program (MIV). The most recent version of tho prototypo data-reductlon program is 
called the Basic Evaluation Program. 

Although tho primary task of Problem 126 has been the development of high 
speed automatic data reduction techniques Incorporating human decision making and 
modifications, additions to the MIV System have consistently been made more elaborate 
than required by the data reduction problem Itself. The results of this seeming 
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extravagance are that a set of elaborste general purpose utility routines are now 
available which constitute the central core of a human-computer problem solving system 
of wide applicability. The development of this system, called SLURP (Servo. Lab. 
Utility Routine Program), la now receiving a major portion of the programming effort. 

SLUR!' Includes all of the facilities of the Comprehensive System developed 
for the Whirlwind Computer by the Digital Computer Laboratory *ta{f, plus a large num¬ 
ber of general purpose routines for extracting and analyzing Information about arbl 
trary computer programs. The governing philosophy of the system Is to allow the pro¬ 
grammer to work on one small part of a very large and complicated problem, virtually 
Independent of the computer which la being used and those aspects of the over-all 
problem which are not of immediate concern. The major features of the system which 
make this possible are as follows. 

1) A group control program which allows the automatic Incorporation of new 
sections of programming Into the system 

2) A manual Intervention (MIV) system of programs which allow the programmer 
to Interrogate and instruct the computer with respect to the over-all problem In terms 
of a specially designed and easy-to-use language represented by cwltches, buttons, 
lights, and visual displays. 

3) A logging program which record* all the MIV actions taken by the pro 

grammer. 

4) An editing program which edits tho logged Information Into an easy-to- 
road record of the manual actions taken, and generates a log playback tap*. 

5) A log playback program which automatically simulates manual action* In 
response to the Instructions on the log playback tape. 

6) An elaborate set of plotting and tabulating programs for data presenta¬ 
tion and record keeping. 

7) A mistake diagnosis routine (MDK) which may be used to abstract arbi¬ 
trary intermediate computed quantities from any program and present these quantities 
In any of the above forma. 

8) Tho SLURP program propor 1* a simulated, generalized, special 
purpose computer which allows the Incorporation of all of these routines Into a 
smoothly functioning system, along with the capacity for continued expansion of these 
facilities. 

All of tho above facilities are under both manual Intervention control and 
programmed control, through the use of appropriately designed special purpose languages 
It Is Intended that the system bo used In the following manner: If a procedure can be 
thought out In advance, then that procedure should be transcribed In the appropriate 
languages so that high speed and reliable execution can be Insured. The automatic 
execution, however, may be overridden at any time by manual Intervention so that new 
Information or new Ideas may be tried out Immediately. Since manual actions may be 
logged, It la possible for a sequence of successful manual actions to be automatically 
Incorporated Into the programmed procedure. By using the variety of data presentation 
routines It Is possible to pursue a line of thought with respect to the problem and 
obtain detailed Information for further study. 
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described hs. 

s=5 rs.sr ^\^^rJs ^ 

will allow accelerated expansion of the system, so that Its p 
be predicted. 

D. T. Hose 

Servoaechanlsas Laboratory 


141 sue SUBROUTINE STUDY 


During the past quarter three programs have been added to the subroutine 


library. 


1 . L8R Fu8a Arc-Sine 

This routine computes the src-slne of an angle In the range -1 4 **1 using the 
following approxlaatlon: 

sin' 1 * - -2 - JTT l a s 1 . 

* i«0 

The routine does not require tho use of buffer registers and cycle counters 
as LSR Fu6 does. Tlie prograa has been tested for accuracy and the error has been 
fouud to be as great as 0.6 x 10' 6 for sin' 1 * very near 90°. 

2 . LSR Fu9 Arc-Slne-Coslne 

This prograa uses the sane approxlaatlon a* fu6a with the added relation 
cos *1' • ^ - sin *x. 

The error can be as great as 0.7 z 10 6 for sin 'x and cos 'x near 90° 

3. LSR MA9 CM Matrix Transpose 

This routine will transpose a square aatrlx which has been stored In core 
neaory. The transposed aatrlx Is left In the sane registers occupied by the 
original aatrlx. 

M. E. Callaghan 

Digital Computer Laboratory 


162 N. DETERMINATION OF PHASE SHIfTS FROM EXPERIMENTAL CROSS-SECTIONS 


Difficulties were experienced In attempting to alnlalze the fitting errors 
In the resonance region. Present and future work will Involve attempting to alnlalze 
the errors of the function In a different area using the saae combination of phase 
shifts previously treated. 
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The prograjaera working on this problem are E. Campbell and E. Mack. 

E. Mack 

Nuclear Science Laboratory 


193 L. EIGENVALUE PROBLEM FOR PROPAGATION OF ELECTROMAGNETIC WAVES 

Hie eigenvalues and eigenfunctions computed previously for the bilinear 
aodel at 50 Me with antenna heights 500' and 30' were used to calculate mode auas at 
17 additional distances. 

Mode sun calculations at 410 Me and 3000 Me using the bilinear aodel have 
been extended to higher mode numbers. At 410 Me convergence has been obtained at five 
distances. These calculations will be continued until convergence at Intermediate 
distances Is obtained. 

Eigenvalues, eigenfunctions and mode sums for an Inverse-square model at 
50 Me have also been computed. Work on these will continue until sufficiently accurate 
eigenfunctions are obtained so that convergent mode sums Bay be calculated. 

The programmers working on this problem are H. B. Dwight and R. M. Ring. 

R. M. Ring 
Lincoln Laboratory 


199 N. LAMINAR BOUNDARY LAYER OF A STEADY. COMPRESSIBLE FLOW IN THE ENTRANCE 
REGION OF A TUBE 

Solutions of the third sat of the differential equations sere obtained for 
the case where temperature dependence of the viscosity, and thermal conductivity of 
the compressible fluid is taken into consideration. 

More solutions are to be obtained for different entrance Mach numbers and 
thermal conditions at the tube wall 

T. Y. Toong 
Mechanical Engineering 


COMPARISON OF SIMPLEX AND RELAXATION METHODS IN LINEAR PROGRAMMING 


In the past quarter several different methods of selecting the zero X , 
which should be made positive at each step In the solution of the transportation' 1 
problem (see Suamary Report No. 45), have been tested. 


Table I gives a comparison of the speed of three different methods then used 
on Whirlwind I to solve a problem with 30 supplies and 260 demands. Methods of selec¬ 
tion require processing of the cost dsta. The first scheme Is the one presented by 
Dantzlg In his exposition of the stepping stone method. It involves searching the 
entire cost matrix to find the zero X which will give the beat rate of decrease of 
the total cost ss it Is made positive. The second aethod, which is the one originally 
coded into the Whirlwind prograa, selects the first zero X which will give a rate of 
decrease of total cost as it is Bade positive. The last technique is a compromise 
between the first two in which the best Xyj for a given plant Is selected, and In the 
following step the best for the next plant is selected, and so forth. 
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1 . 

2 . 

3. 


Table 1 Co.parl.on of Method. - 30 * 260 Transportation Probl, 


Method 

Best X 

I--- 

Number of 
Iterations 

No. of Examinations 
of Entire Cost Matrix 

Total 

Time, Min 

508 

508 

19.4 

U 

Firs' X 

2200 

27 

25.0 

U 

Best E tJ for 

a plant 674 

44 

9.6 


It la evident that the co.pro.lae aethod la better by 
either of the other, and, hence, .111 be used for future rune, 
of the routine la contemplated. 


a factor of 2 than 
No further development 


3. B. Dennla 
Electrical Engineering 


225 B,N. NEUTRON DEUTERON SCATTERING 

The calculation of n-d acatterlng length. with an expanded 4-parameter trial 
function haa been performed for the apeclal caae M =0 (no polarization). The result, 
differ very little fro. those of the earlier two parameter trial function. The next 
step .111 be to u.e non-zero values of the polarization para»otor. 

The result, of the calculation to date are being presented at the New York 
■eetlng ol tho American Physical Society.* The programmer Is H. Paul. 

• Bull. A., Phys. Soc. II 2, 71, 1937. 

L. Sartorl 
Physics Department 


231 B, N. AN APPLICATION OP MONTE CARLO METHODS TO NEUTRON DIFFUSION 

Progress .as made In trouble shooting the assembly of subroutines, but some 
errors still remain. 

M. Troost 

Chemical Engineering 


244 C. DATA REDUCTION FOR X-I FIRE CONTROL 

Through the use of new data correction achieved by modification of a program 
tape previously prepared by Dr. J. If. Stark of the M.I.T. Instrumentation Laboratory, 
data reduction on this problem came to completion during this quarter. The data has 
nos bean processed and Is In final form. 

No more computing time is needed on this problem; however, the program 
tapes are being saved for possible future use. 


S. J. Madden 

Instrumentation Laboratory 
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245 N. THEORY OF NEUTRON REACTIONS 

More test runs .ere made In tho vicinity of the P and D waves In order to 
locate a t and c that will give a more satisfactory fit to the experimental values than 
5 = ,08, c a 1.65. When such values are found a complete case will be done. 

E. Campbell 
Physics Department 


253 N. AN AUGMENTED PLANE WAVE METHOD FOR IRON 

In the augmented plane wave method of handling tho one-electron problem In 
crystals, It Is necessary to Integrate tho radial wave equation In the chosen spherics) 
potential. Tho accuracy of the numerical method (Noumerov method) Itself and the ac¬ 
curacy of the Whirlwind routine for carrying out this method have boon under Investiga¬ 
tion. Use has been made In these checks of hydrogen-llke potentials where the solu¬ 
tions are analytically known In the form of Coulomb wave functions. 

At present It appears that the accuracy of the extend Integration routine la 
quite sufficient for our purposes. 

The other routines for the APW method (developed by M. M. Saffron) have been 
tested and It appears that one should be able to carry through a one-electron energy 
band calculation via this method In the very near future with an accuracy limited prl 
marlly by the approximations of tho one-election method Itself. 

J. H. Wood 

Solid State and Molecular 

Theory Group 


257 C. HORIZONTAL STABILIZER ANALYSIS 

Some production rune havo boon made with the master program. There are a 
few minor uncertainties In the problem under Investigation which may dolay the continu¬ 
ation of the production runs. Once these uncertainties have been eliminated, the pro¬ 
duction runs shall bo carried to completion. 

The programmers working on this problem are N. P. Hobbs and E. S. Crlscione. 

E. 8. Crlscione 
Aeroelastlc and Structures 
Research Laboratory 


261 C. FOURIER SYNTHESIS FOR CRYSTAL STRUCTURES 

In the past quarter, the two-dimensional Fourier synthesis programs have 
been used to develop a theory of partial Fourier synthesis. In this kind of synthesis 
only part of the total number of Fourier coefficients are used. The resulting synthesis 
Is useful when dealing with crystals which can be said to have a substructure. For 
these crystals, the "ourier coefficients can be split Into two sets, not necessarily 
equal In number, one of which contributes to the electron-density pattern of the sub¬ 
structure, the other to the electron-density psttern of the complement structure. 
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Since the complement structure Is ordinarily swamped out by the substructure, especially 
In Patterson syntheses, this method has proved valuable In investigating crystals 
ehlcb can be referred In any way to substructure theory. 

M, J. Buerger 
Geology Departaeut 

262 N. EVALUATION OP TWO-CENTER MOLECULAR INTEGRALS 

Soae Couloab Integrals wore evaluated by Dr. P. J. Corbato's program. 

The results were compared with those computed previously by an alternative method. 

The results differed, In some cases agreeing to si* decimal places and In othera two. 

H. A. Aghajanian 

Solid State and Molecular 

Theory Group 

264 C. OPTIMIZATION OP ALTERNATOR REGULATING SYSTEM 

The method for minimizing a nonlinear function In the presence of nonlinear 
equality constraints by the methoo of steepest descent has been extended to the In¬ 
equality constraint problem. This technique may be used to solve the problem of de¬ 
signing an optimum system, when said problem Is phrased as a nonlinear programming 
problem. Trial computer runs have been made using the simplified aircraft alternator 
representation described In the last quarterly report, and have shown that certain un¬ 
desirable characteristics of the method are emphasized when Inequality constraints are 
considered. 

In order to correct these faults In the program, the analytical study of 
the techniques used therein has been Intensified. This study should be sufficiently 
advanced by the next quarter ao that coding for an appropriately aodlfled program nay 
be Initiated. 

Programmers working on this problem are J. B. Dennis and R. R. Brown. 

R. R. Brown 

Energy Conversion Laboratory 

273 N. ANALYSIS OP AIR SHOWER DATA 

There have been minor changes In the program for the purpose of finding 
systematic errors In the air shower data. The error analysis Is aided considerably 
by the fact that the parameters of Interest are highly ov.rd.termlned by the set ^f 

th * * ,r ,ho " r — 

reduction. ^About ^Doo^showers t *ave°been*processed"and ^rkT .T tlTt^on 

tinu. during the next quarter. Progra-er. working on t III ^bl- a^lc^d 


P. Scherb 
Physics Department 
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» »» .m, „ aPmiou . 

- •<-• p W ... 

-...............*“ 


M. Karakashlan 

Nuclear Science Laboratory 


ENERGY LEVELS OP DIATOMIC HYDRIDES L1H 


hydride molecule using the electronic wav. f^cUon^hl h* Pr0p • ^tl • , of ,b « lithium 
the past quarter we have calculated the dipole moment !* h ‘” e obt «l"«<l. During 

nuclear distance, for which we have f^ctl^T * ,UnCU ° n ° f lh# lnt " 

which decrease, to zero Infinite separating ' , * F * • h0 ''* * “‘P 01 * BO * ent 

a maximum , t an Int.rnuclear separation .bo^one .?Lr r ° •• p * r,tlon **>‘ch reaches 

briu- .. A1 taough the c^cullted ^ueTs th “ thB «•-»»»- 

mental value, the shape of «he curve (dipole m«entTr ' h,n ,h ” 

pos.rion of the maximum ... possibly qul£ Ic ^" ‘■t.rnuelm.r distance) and the 
verification. ^ accurate. This await, experimental 

rating llthlul” u T‘\T a,U ’. orblt * 1 treatment Incorpo 

binding energy and other’-oleLui'p^ti::^ £ ‘w ln th « 

angular correlation. operties will result from some Inclusion of 

Programmer. working on this problem sr. A. M. Karo and A. R. Olson. 

A. M. Ksro 

Solid State and Molecular 
Theory Group 

288 N. ATOMIC WAVE FUNCTIONS 

completed. Extensive Port'll.™ been^ono^irh^r^M * d ^ “ rlUp r ‘ ,port * has been 
(see January 15, 1957 Quarterly Progress Kenort f fUn »i l0 c* ln v * ^lou • environments 
Theory Group,. Th.s work h.s^^t'X^" '"l ^ •" d 

spent on Cl function, using the same number “° Unt of ‘ h *« b »®" 

stabilizing environment We intend r P *** ter “ ** ln the hut without a 

tlons: 4S2 3d«, 4 S 3d7«nd 3 dS „ lth th^o“«t“lo U ^“o n, *‘2 thre ° confl « ur * 

- -- itr^r^r io i -"'“"-'-.m 

an^uivaLnt calculation and the limited amount of Itor^n 


Programmers working on this problem sr. R. K. Ne.bet and R. E. Wst. on . 


R. E. Watson 

Solid State and Molecular 
Theory Group 
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290 N. POLARIZABILITY EFFECTS IN ATOMS AND MOLECULES 

The mathematical and computational aapecta are the same as previously given. 
A summary of our results to date and an analysis of the physics of the problem as It 
relates to electronic save functions for the hydrogen fluoride molecule Is given In 
the January IS, 1957 Quarterly Progress Report of the Solid State and Molecular Theory 
Group. 


L. C. Allen 

Solid State and Molecular 
Theory Group 
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rt. sr- t jssss? 


Ihla work la sponsored by the 0 s 
Programmers working on this problem are J. S. 


Air Force under Contract AF 33(616)2392 
Prl Ege. Jr. and L. K. Berman 


J *. Prlgge, Jr. 
Aeropbyslcs Research Group 


312 L. ERROR ANALYSIS 


300 L. TROPOSPHERIC PROPAGATION 

The present phase of the Fresnel integral computations have been completed. 
The autocorrelation function and power spectrum of a sample of 2290 Me. field strength 
data have been determined. A subroutine for generating Incomplete elliptic Integrals 
of the first and second kinds has been developed. 

An extensive Investigation of statistical methods aultable for analyzing 
data obtained In the Aurora research program will be carried out. 

Programmers working on this problem are W. C. Mason and P. A. Duffy. 

W. C. Mason 
Lincoln Laboratory 


306 D. SPECTRAL ANALYSIS OF ATMOSPHERIC DATA 

During the quarter most of the computations were completed. It la antici¬ 
pated that the problem will be terminated during the next quarter. 

B. Saltzman 

Meteorology Department 


309 B.N. 


PURE AND IMPURE POTASSIUM CHLORIDE CRYSTAL 


to time. 


All programming has been completed. 
Programmers working on this problem 


Production runs are made from time 
are E. Hutcheson and L. Peterson 


E. Hutcheson 
Lincoln Laboratory 


r 


It la planned to refine the statistical i 
derivatives from suitable responses by all methods 


lethod and to try to extract lateral 
being Investigated. 


Mazzola, 


Programmers working on this problem are T. M 
M. N Springer and L E. Wilkie 


Carney, M 


Hoult, L L. 


M N. Springer 
Aeropbyslcs Research Group 


During the past quarter, the results of the problem were prepared for publl 
. l0n '»K TO ° T ’ obt,lned suggested several additional details to be investigated 
During the coming quarter, calculations concerning the charge density of the perfect 
crystal and an Investigation of the conduction band are planned. ** 


327 L. PREDICTION ANALYSIS 


The Maximum Likelihood 
No. 46, has been programmed. It 
next few weeks. 


Criterion Mothod, first described 
la hoped that checking out will be 


In Summary Report 
completed In the 


L. P. Howland 

Solid State and Molecular 

Theory Group 


„ Pr ogrameere working on this problem are E. Hutcheson. 

* • «iii uno 


L. Peterson and 


310 C. TRAJECTORY CALCULATIONS FOR A ROCKET DURING POWERED FLIGHT 


In the last quarter, the program was used 
parametric study of powered-flight trajectories of 
this study are presented In the following report: 


extensively In carrying out a 
long-range rocketa. The results of 
Prlgge, J. S., "A Digital Computer 


E. Hutcheson 
Lincoln Laboratory 
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329 N. FIRST APPROXIMATION SOLUTION ON ORE BODY 

The problem of determining the optimum ilie and shape of the fundamental 
block to be used In the Inverse Interpretation of field resistivity data has been 
solved and production runs of the various geometric configurations have been made. For 
blocks at or very near the surface of the earth the solution Is poor but because of the 
fact that geologically a layer of overburden of resistivity comparable to the back¬ 
ground value always exists, It Is unnecessary to use such blocks In the Inverse prob¬ 
lem. 

Tlia cross-section of a general resistivity profile for pole-sender and 
dipole-receivers has been established and attempts have been made to Invert the re¬ 
sulting matrix of conductivity coefficients and geometric factors. Thus far, there 
has been little success primarily because of the fact that In the use of the Crout 
Method a great loss of accuracy occurs If the diagonal elements are not greater than 
the elements to their right. Usually It Is possible to rectify this situation by re¬ 
numbering the equations, but this has not resulted In sufficiently Increased accuracy 
in this problom. It will probably be necessary to determine the corrections to be 
made from the matrix formed from the residuals, and It Is hoped that convergence of 
this Iterative procedure will result. The correct Inversion of this matrix will 
enable a rather detailed Inspection of the distribution of the conductivities of the 
earth useful in mining evaluation. 


N. F. Ness 

Geology and Geophysics 


334 C. PARAMETRIC STUDY OF COUPLING AND DAMPING 

IXirlng this quarter, a few production runs wore made to Investigate the 
possibilities of using this program to Investigate airplane flutter. 

K. R. Wetmore 

Aeroelastlc and Structures 
Research Laboratory 


337 N. NONLINEAR 2ND ORDER DIFFERENTIAL EQUATIONS IN THE THEORY OF ELASTIC SHELLS 

During this quarter the main body of the program has been rewritten It 
was necessary to split the program Into parts, using auxiliary drum group, in order 

;rr: s.sir.’rss.ss irs.'u-; »:s ,u ; - *“■*- 

1 lmlnary computation, are desirable, which are -xpected to give ° ‘ \ ! Pr * ‘ 

make the 704 program to be written more eff.ctlv.lv Vl , , tC 

carried out In the near future. Computations on the «,lr”lnd ^ 1,111 ^ 

pectod to be completed by the end of February. rlwlnd I computer are eg- 
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H. Welnltschke 
Department of Mathematics 
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341 C. 


STATISTICAL AND DYNAMIC FORECASTING METHODS 


report, were^uccesafulljT.“n^ to dilii^b? I*"!!" pr#dlct,on d «crlbed In previous 
Phere. The result, of the “ ‘a* 5 °° * b ‘” Vel ,n th * 

the coeff leient.^'of”the^emp'rleal orth^oMl'f "" m* P r ««ure map based on 

on (1) observations of sea-level pressure on two'. 100 * conUnued - Prediction was based 
observations of .e.-i.v.Tpr." e ir. pr ::dTe::u^:: d c ::!: i r: «?*-—’ 

coefficient, a. predictor, th. first “«u!n2d5? * d * y ' " lth 

cent of th. variance of ...-level pre.sur^t^ fon^wlel T C "‘ * n<1 ^ • econd 63 

sampled* 6,1 ing c^ «'o d t«?: t ; 0 h ^*,„ b ::^ n o d b : n * t tne<i ^ 




E. Kelley 

Meteorology Department 


343 C. WEATHER PREDICTION 

Austin, Dep.^eT ^iT ” r * " ltten “ nd - r th# ■ up "”‘— d 


E. Kelley 

Meteorology Department 


346 B. COMPLEX SPECTRUM ANALYSIS 

s-s-mZ ^,r,. s “sr 
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This urlci thi first time thit work of this sort h»l been cirrlid out on 
digital computers. Thw computer hu proved to be in extremely valuable tool In the 
analysis of spectra, and It will probably be used extensively for this work In the 
future. 


J. W. Lindner 
Spectroscopy Laboratory 


350 D. LONG RANGE FORECASTING 

The computation of eigenvalues wentloned In the last report has been com¬ 
pleted. Lag correlations for temperature and pressure forecasting have been obtained 
from them, as have contemporary correlations for temperature specification and pressure 
map reconstruction. To complete the problem, pressure forecasts still have to be 
computed. 


Prog rummers working on the problem are B. Shorr and D. L. Gllun. 

D. L. Gllun 
Meteorology Department 


334 D. RESPONSE OP A SINGLE STORY CONCRETE BUILDING 

The effect of various paraaeters Is being Investigated. It Is expected that 
the problem will be completed early In 1937. 


S. Namyet 

Civil and Sanitary Engineering 


361 B,N. GROWTH OF FATIGUE CRACKS 

The integral equation described In the last progress report has been changed 
In order to determine the rate of crack growth from notches of finite radius. This 
new equation has been programmed and Is currently being debugged. Upon completion of 
the debugging, about 5-10 runs of about 3 minutes each will be necessary to cover all 
parameter variation. 


The finite difference solution for another plasticity problem (the torsion 
an axially symmetric notched bar) la now under way, using Whirlwind. A solution for 
one set of parameter values Is being done unually, using relaxation techniques It 
Is planned that the Wnlrlwlnd work will be completed by Hay. 


J. B. Walsh 

Hechanlcal Engineering 


364 C. BLAST RESPONSE OF ROTOR BLADES 


Routines for calculating the approxlmat 
a finite set of assumed nodes have been completed 

In the future, these "normal modes" wll 


norul modes of rotor blades from 
be used to Integrate the equations 
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ot motion of the rotor blade under various blast loading.. 

A. Holland 

Aeroelaatlc and Structures 
Research Laboratory 


the transient h«rtrL*f^r 0 problel , '.fr i p,!!^- ‘q <l f^*fwsntlsl equation, describing 
Magnet lc Drum Calculator and the Whirlwind Computer « !J”" 0n b ° th lhe 650 
Quarterly Report. During -he present Ascribed In the preceding 

for the entire system u.!n« s.tetn dl»t™t, ' l "T” ° f •° lutlon> •<”•« obt.m.S 
value. * ■ ever#1 different physical confIgur.tlon. and parameter 

and these. In turn, serttd tf^httkTht °" J h * 6S ° ‘ L, * n ' tlc Bru " Calculator 

ceeding case studies on the Whirlwind Computet" Tl tn * r f* Ult * ob, * ln * d °n *uc- 
-gree to 7 and In aome ca.ea 8 detl«, TiluY/ ^ r0 ' UU ‘ ,h ° computer. 

-unwind achieve oT^i^^T^ '“iTtV 

the characteristics of this compiler. Th ® following section describes 


Special Programmed-Arithmetic Co 


apller for Whir lwind 


Heat Tr.n.^°pr“L r T^ ... - 

to a maximum. J tl e of Increasing the speed of solution 

The ba.lc compiler 1. composed of th. following routines: 

ulator by III ,‘" U th * . .*««.- 

£.2ir ult ,n th * u^^^T-Str- 

ar .uhout'‘JistssHh; ^i:. h :in c , OB t t h e e nt :..°u , it lh :n‘:h c r utor , to th# «* 

“ 

exponents differ by ■ ore than *1 n^nor-eT 1 * 10 " °' ^ ““‘h”" ' h °‘ # 

The folle.tnr3^ r i5 f rtuU^rtive r0 ^"ru"n b, ' en ^ b ‘" C C " pll#r - 

..... „is -sr-srs! 

.........' rr. r. -jr z,rr.r„r... 
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(c) Printout - TNI* routine print* (without disturbing the content of the 
accumulator) on Tape Unit No. ] or I using 8 decimal digit mantissa and 2 decimal 
digit exponont. The print routine used the multiplication routine 

The Scope coordinate routine plots the divided coordinates on the scope. 

It does not use the compiler. Point plotting routine unpackages the number in the 
accumulator and plots It In correct position on the scope. The print routine does 
not disturb the accumulator or use the compiler. 

All the routines operate on the accumulator and the AR therefore the de¬ 
sired number* should be placed there before transfer of control Is effected. All 
routines operate on normalized (scale factored) numbors. The Insertion of an un- 
normal lzed number Into tho compiler will cause a divide alarm. 

In order to save time, there are no bullt-tn checks In the routines, and 
all possible alarm* are regular Whirlwind alarms. The following alarms are possible: 

(a) Arith metic ove rflow If the binary exponent of a number gets larger 
than 2 15 or smaller than 2~ 15 . 

(b) Divide alarm m case an attempt Is made to divide a non-zero number by 
zero. (A zero divided by zero will result In a zero). 

Tho routines are written to provide the fastest possible speed and no at¬ 
tempt was made to minimize the number of instructions. 

The comparative speed and length of the basic routines are given beloe: 
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Length of 
New Routine 


Multiplication 

Addition 

Subtraction 

Division 

Clear and Add 

Clear and Subtract 

Store 

Jump 

Conditional Jump 


0,645 - 0.500 
1.000 - 0.200 
1.000 - 0.200 
1.217 - 1.150 
0.111 
0.111 
0.111 
0.015 
0.030 


In order to Increase the speed of operation, the major and minor parts of 

J7 b T,h he • CCU “ UUt ° r tnd ,h0 AR been placed In the flip flops. 

The availability of the above compiler reduced the machine solution time bv a factor 

» prevlously ,rt * ln C “°‘ " ' hOTO PTfo»ed 8 time. a. till 

Commit. , P "* r T“ r ‘ W ° rklng ° n thl « P robl « "" *. Zabludowsky of the Digital 


DRSIGN OP SPHERICAL SHELL SEGMENTS 


Stresses In a spherical shell segment due to7t7‘ ""' ,nU,n > sieved •"» »U 

the future one program only must be slightly 57 h ‘ V * b '* n >» 

the numerical result, for the present problem wolutlon will complete 


K. Traum 

Civil Engineering 


COVERAGE ANALYSIS 


Run. are no. being mlde!‘ n ?h'7p7grM I7o£lTtent«t J,7"h Pr °* r, “*“ 1 ‘ nd ch8cked °ut 
be extended or modified! °" Y "owever, a. lh« Analysis nay 

B. Hutcheson 
Lincoln Laboratory 

380 B. SWITCHING CIRCUITS 

switching functions of^is or’Je^r v^Ubf«'Si B S.jr , c«pl"t P d OC *^ r, f f ° r ,1 " pUfyln « 

==ilsl^HHHH:==r£l 

C. Roth, Jr, 

Electrical Engineering 

382 B. CA ISOLATION OP PRIME NUMBERS 

g.aa I h.« C no! iZ 777* ° U * , “ P "’ copl8 * of the pr ‘« calculation pro- 

77: a completed my work on the Whirlwind computer. In .bout three week. 

I intend to eub.l, . report telling exactly whet e.ch of my progr... doe. P„7!7 

7* T ° il , 77 dB t ° d ° ,U 7 her work on Protr^m., 1 suggest reference to my the.li 
advice evelleble during the semester, Pebniary - June, 1957, for perion.i ' 


P. M. Verzuh 

Office of Statistical Services 
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k-t! .nrrsr; =s “.r “ LI—r— 
ssrs ia-* «»— 

In my thesis. 


H. I. Cohen 

electrical Engineering 


383 C. STOKES PARTICLE VELOCITIES 

Experimentally obtained parameters were fed to WWl and all data baa been re 
duced and plotted, Further work will conalat of running only a few specialised para¬ 
meters through the computer to aid In the writing of a final report. 

Programmers working on this problem are T. A. Harlow and J. 0. Housley. 

T. A. Marlow 
Hydrodynamics Laboratory 


385 B. FEED PLATE LOCATION 

The range of optimum feed plate locations was determined for each of 
twelve ternary distillation systems at five different values of the reflux ratio for 
each syotem. The data obtained from these computer calculations were sufficient to 
permit an analysis of the factors Influencing the optimum feed plate location In 
multicomponent distillation systems. 

Pending final acceptance of the thesis written using the data mentioned 
above, the problem Is completed. 


I. Rluard 

Chemical Engineering 


386 C. FREE CONVECTION PAST A VERTICAL PLATE 
Problem Description 

The total differential equations (representing momentum and energy conserva¬ 
tion) to be solved simultaneously are the following: 

Note that O' refers to differentiation with respect to the Independent 
variable, »| . 

( t ••• + (1 +• q) ff" - (I + 2q) r a a - Bh 
) (Pr' 1 ) h” + (1 + q) fh’ - (1 + 4q) f • h . 0 
subject to the boundary conditions: 
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* °* f . f . 0, h . 1 

1\ r- h • 0 

It is proposed to employ the 0111-Kutta routine for Integration purposes 
**““*•• wlU J** “ d0 of f (0> * nd h <°> ln Umu °f f W) and h(«e) and the integration 
will be repeated until satisfactory and point conditions are met An Influence coeffi¬ 
cient routine will be omployed to enable the computer to determine a successive se 
quence of f (O) and h (0) values which will converge upon the proper values on the 
basis of successive f (“) and h«e) results. This technique was used previously by the 
present programmer for two-point boundary condition problems. 

The scope of the computing problem Is defined by the following: 


(1) 

0 s >1 i 5.0, l.e., wj ■ 5.0 

O oo 


(11) 

» 0.1 for Integration 

Increments 


(111) 

q * -.10, -.05, or -.01 ) 



(lv) 

Pr > 0.01, 1.0, or 100 

present Interest Is ln 11 

coablna- 

(v) 

B > 0.01, 1.0, or 100 

tlons of the constants q, 

Pr, and 

(vl) 

Results to be printout (presumably on magnotlc tape) 
of f, f, f , h, h' as functions of for each case. 



The required estimates of f (0) and h'(0) should be fairly accurate since 
exact aolutlons are available for q - -.25 and -.20 from which the approximate values 
needed here may be extrapolated for the q = -.10 case, etc. 

Progress 

Several similarity solutions tor the dynamic and thermal boundary layer 
profiles have been obtained during the past quarter. The solutions correspond to a 
family of surface temperature distributions of the form 1 + Bx <J whero B and are 
constants and x represents distance along the plate Several additional values of B 
are contemplated for future computatlona. All results are presently based upon a 
Prandtl number of unity. 

J. R. Baron 

Naval Supersonic Laboratory 


387 C. DETERMINATION OF VELOCITY POTENTIAL 
Problem Desc rlp t Ion 

In order to evaluate the unsteady lifts and moments acting on wings under 
various loading conditions at supersonic speeds, the wing planform is divided Into a 
grid of elementary areas, or boxes. The velocity potential Influence coefficients 
(VPIC) evaluate the velocity potential at each given box due to unit downwash at another 
box. The VPIC's are then properly weighted and multiplied by the actual downwashea to 
obtain the velocity potential distribution on the planform. The velocity potential 
Is suitably Integrated to determine the generalised forces acting on the planform for 
each Mach number and reduced frequency. 
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ft. portion of tb. problem th.t Mil ««»»*• ** ^ 

in, matrix operation., .her. most of th. numb.r. ln»olv.<l ar. complex. 


I. 


Weighting of th. oo.fflcl.nt -trie. t.rm by t.rm, to account for 
th. .ffect of partial box.a on th. planform 


C DD : C DP 


r dd , r dp 


S DD ; C DP 


X 

- - - - — 

:S 

• _ 

Sd ; c pp 


, r pd 1 *PP 


l C PD 1 PP. 


A. 


3. 


Solving ..qu.ntl.Hy for th. down.a.h.. In th. alxed region. n.ar th. 
edge. of tho planform: 

{M * * [ C Dd] 1 [ C D(j{ W p} 
wh.r. [C DD ] 1 . a diagonal matrix. 

Solving for th. velocity potential dlatrlbutlon on the planform 


■[ C roi c Pij^ p | 


.here: D - mixed or diaphragm region 

P = planform region 

C • unweighted velocity potential coefficient. 

C . weighted 

R = weighting factor, to account for the effect of partial 
boxe. 

W = downwa.h 

if a velocity potential 

A total of about 18 ca.e. will be investigated. The large.t matrix will 
be 46 x 73 complex. 

Program. During th. Pm.t »iarter 

A generalized routine ha* been formulated for complex matrix multiplication 
capable of handling aatrlce. up to order about (18 x 18) for the premultlpll.r and 
(18 x 1) for the po.tmultlpll.r (or equivalent) where all the number, are complex. 

A generalized routine ha. been formulated for a aequentlal solution of com¬ 
plex equation., c.pablo of solving up to about 18 equations. The matrix of coeffi¬ 
cients 1. triangular and all the number, are complex. 

Por determination of supersonic velocity potentials on wing planform. a 
generalized routine ha. been devised which consists of 2 long tapes: 

a) Input tape which store, the basic aerodynamic Influence coefficients, 
th.lr appropriate weighting factors, and a given mode shape on th. drum; 

b) Computation tape which take, the input data and compute, the resultant 
velocity potentials corresponding to the given node shape. 
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“V* to yl.ld satisfactory velocity potential. corr..pondtng to two mod. shape, of 
delta planform at a given Mach Number and reduced frequency of oscillation. 


Future Plans 


Additional computations will be performed utilizing th* existing routines 
for determination of velocity potentials for the delta wing, corresponding to many 
more combinations of mode shape. Mach Number and reduced frequency. 

The existing routines will be modified (the aodlflcatlon. will be relatively 
minor) to handle the velocity potentials for a straight wing In supersonic flow and 
corresponding calculations will be performed. 

fte program will be modified (these modifications are also expected to be 
relatively simple) to handle swept wings and corresponding calculations will be 
performed. 


S. Gravltz 
R. Stapleford 
Aero.la.tlc and Structures 
Research Laboratory 


388 D. TEMPERA TURK DISTRIBUTION IN AIRCRAFT GENERATORS 
Problem Description 

Present aircraft generators sre cooled by directing a blast of air through 
Internal ducting. No analytical Investigation of Internal temporature distribution 
bas yet been aade, with the result that placement of ducts Is somewhat arbitrary, and 
prediction of machine temperature capabilities Is Impossible. 

Th* present problem represents an attempt to utilize relaxation aethods for 
determining generator temperature distributions, along with numerical Integration of 
fluid flow equation, to calculate.the heat-transfer coefficient In each part of a 
machine. For solution of the entire problem, It Is planned to use the superior speed 
and storage capacity of the 704, but Whirlwind will be used for preliminary ord.r-of 
magnitude" fluid flow calculations, intended to determine which flow obstructions are 
most significant. 

Theory may be idealized as a calculation of the flow In a tube, with heat 
-ddltlon and friction. The appropriate equations are Integrated by taking small 
enough step, so that an average fluid temperature may be used. 

Progress 

During the past quarter there has been substantially no progress because 
of difficulties with the underlying theory. 


R. Moroney 

Servomechanism. Laboratory 
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391 L. MAGNETIC RELAXATION IN THIN FILMS 

Problem Description 

Hiln Mine (- lOOO A) the 

bistable ■•■ory el •■•ms for digital computers f , which is sub- 

relaxation tin. -1th coincident current 1. of the order of 10 £ 

.tantlally fa.t.r than present memory eLeenta. Kxperl^n. to. .^ 

"S”r!^^ , c«'i.nd 1 ir“no“ d .^quI»"i«r*y baaed on the Lmnd.u-Uf.hlt, equation 
of"aotlon*of > 'th^oagnetlxatlon can only be treated by nuaerlcal method. on a -chin, 
such as Whirlwind I. 

Th. specific problea 1. to solve the pair of couple differential equations 
8 = eg Y + r(0) 

Y= -ay + g r( 0 ), 

with r(0) - h sine 4- hj cos e sin 9 sin 8. and a and g known constants of 
the magnetic material. 

Solutions aunt be obtained for a set of points In "h-space", 0 ‘ h g ‘ 2 and 
0 <hi<l. It Is proposed to use the program of Lanlng and Zlerler which translate 
equations directly into Whirlwind I language. 

Progress During the Past Quart er 

Lanlng's prograa was used and was found to be too slow. The problea has 
been reprograaaed using the Comprehensive System. 

Programmers working on this problem are J. Frankovlch and D. O. Smith. 

D. 0. Smith 
Lincoln Laboratory 

393 N. THE INVERSE BREMSTRAHLUNG SPECTRUM 

The yield Y(E ) of a photonuclear reaction produced by a bremstrahlung 
spectrum S(B ,k) Is meaSured as a function of the maximum energy E of the spectrum. 

It is usually desired to know the cross-section c(k) tor the reaction as a function 
of photon energy k. 


Y(B_) 


I* 

f 


S(B ,k) c(k) dk , or approximately 


Ik* 1/2 


YlQet- 1/3)0 - t 


z 

h=l 


S< [n + 1 / 2 )t, nt) c(nt) and 
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S '(nt, (m + 1/2) t> Y( [a + l/a) t). 

S Is the inverse of the triangular matrix S. 

Using the Schlff bremstrahlung spectrum for 8, Whirlwind will be used to 
compute the elements of S for ten twenty-by-twenty S matrices. Three of these will 
cover th. range 5-25 Mev with t . 0.50 Mev and the other seven .111 cover th. rang. 
8-18 Mev with t * 0.125 Mev. 


c(nt) - t* 1 1/2 

£ 


H. Paul 

Nuclear Science Laboratory 


394 C. AUTOMATIC PROGRAMMING FOR NUMERICALLY CONTROLLED MACHINE TOOLS 

The completion of the MIT Numerically Controlled MlUlnE Machine In 1952 In¬ 
troduced a nee era In the automatic translation of design requirements Into finished 
parts. The primary problem associated with thie type of machine tool, the preparation 
of detailed Instructions for making s part, has been studied using the Whirlwind com 
puter. Under problem 132 a library of subroutines was developed to assist in data pre¬ 
paration for the MIT Numerically Controlled Milling Machine (see Whirlwind Summary 
Reports 33 through 43). Under problea 250 the feasibility of applying sodern automatic 
programming techniques to machine tool programming was successfully demonstrated for 
plane parts bounded by straight line and circular eegaents (see Whirlwind Summary Re 
ports 41 and 42). Problea 394 Is concerned with the extension and elaboration of these 
automatic programming techniques to arbitrary surfaces which can be machined by 3 and 
5 axis numerically controlled milling machines. 

The problem may be viewed as the design of an APT System (Automatically 
Programmed Tool System) combining the human operator, the machine tool Itself, and an 
APT computer, which la a combination of a general purpose computer and existing machine 
tool directors. The objective Is a system of Whirlwind computer programs which will 
simulate an APT computer with a language which Is convenient for human expression of 
machined part designs. A design expressed In this language will be automatically 
translated by the APT computer Into machine tool Instructions. 

Whirlwind programs have been written which allow programming of a part by 
means of arbitrary curves In the surface of that part. At present the system is re¬ 
stricted to planes and spheres but may be arbitrarily expanded. Work le now In pro¬ 
gress on an advanced system which will sllow programming by means of entire regions of 
the part. Following completion of this advanced simulated computer, attention will be 
directed to the determination of the moat appropriate English like coding language for 
communication with these APT computers. 

Associated studies are being made for the determination of appropriate 
languages for the feedback of Information from the APT computer to the human operator. 

A program has been written by which the Whirlwind computer can draw on the output 
oscilloscope an axonometrlc projection (near perspective) of the part being made. The 
program allows arbitrary orientation of axes and arbitrary scaling so that any view of 
any portion of the part can be shown. The program has several modes of operation, but 
at present the picture Is displayed as the path of cut vectors swept out by the motion 
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of the cutting tool over the surface of the part. Engllsh-llke computer to human 
languages "111 alao be devised. 

The programs are currently being developed by Douglas T. Rosa, Harry E. 
Pople, Jerome R. Wenker, George Zames, Samuel M. Matsa. Servomechanisms Laboratory 
staff members. This work Is sponsored by the Air Materiel Command through Project 
6873 In the Division of Sponsored Research. 

D. T. Ross 

Servomechanisms Laboratory 


395 L. PHASE ERROR CALCULATION 


The problem Involved la one of calculating the far-fleld energy distribution 
of a radar antenna that has the form of a parabolic torus. In using an antenna of this 
configuration there Is a certain Inherent phase error problem whose effect can be de¬ 
termined by constructing a model and measuring energy dlstrubutlon patterns. As an 
enhancement of this method, one may also calculate tho energy pattern of the antenna 
An expression for doing this Is as follows: 


r b/ 


2TT) 

s ». (A + 


slnO (x cos 0 + y sin 0 ) 




angle . ( 0 * t 8 Z'VTFl'l 

in Step, of 0 . 01 . HmpJt ti" for .U x in «.ro?VS, rOU, l*!‘ ” ( -°-* 4 5 * = 0 2 

contributions fro« all * and y at each f or . ' 1# Havln K ■uMed the energy 

step, and repeat the above Ch, "‘ a # ^ 01 


Number of Iterations: 

y - 48 
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0 180 


Total number of Iterations: 


48 x 48 , 50 , ,80 ? 20 , ,0 6 Iteration. 


incorporated l^cIlmUtlo2ii*of ’the^dl.Uon" h * Ve h*”" ° bt * ln " 1 The, » »m be 
antenna. ° r th * r,d “ ,lon pattern of a toroidal reflector type 


D- C. MacLellan 
Lincoln Laboratory 


SUBROUTINE STUDY 


to invert * Ut9r ~* «‘ b ° d ' *»<>*" 
stution), Craig'S, Crout's and Frovenlus'' h *" ‘ h * • llBln * ,lon (u«lng back subatl 
m particular, *.e have I 6 Tl ^ 'lU^ondlM —'““.oned -trices, 
would like to do is: complex. Ill conditioned matrix to Invert. What .. 


1) Find the Inverse of the particular matrix: 


conditioned’.a trices" * • Ubr ° UUn ' ‘—rporstln, the best ' method available for 

order 1 . t~ ^g" tl‘t IWnto".^*^ •»» -trice, whose 

some sort of partition!^. ' * “ y P °“ lb,y "*>• *» *> tM. by 

During the past quarter, the problem waa programmed to solve . u . 

J. Kuseaan 

Digital Coaputer Laboratory 


397 N. RESPONSE FUNCTIONS OF AIR SHOWER DETECTOR ARRAYS 

nrehl . ^ ! hOW * r “ ,re to detected by a specified array of detector.. Tho 
proWem 1 . to determine the efficiency with which the array detects air shower, of 

* b °“** e ' f " U wlthln • “«•»■ •>•«. A and -hose arrlva. direction. 11. .man 
* certain solid angle u>, 

We cal1 R<N9)dNde the rate of detection of showers with sizes between N and N 
+ dN and zenith angles between 9 and 6 + dO. Wo call S(N0)dNd«» the rate of showers 
with sizes between N and N + dN which arrive with a zenith angle 6 within the solid 
angle du> and cross unit area perpendicular to the direction of arrival. it I, related 
to S by the equation 
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R(N8) * 2T (In 0 co* 9 S(N9) 


/ dr / to / « CN.«.r.a.f ) 


where 


t * TT Pj P (H.O.r.a. f 


r . 


h . | u - */ + (y - y/ « t > + ■ <» - *i>] *} 


1/2 


i . tin 0 sin f , 


■ • sin 0 co* p , 

p (a ,f) • the probability that the 1 th detector Jflth projected are* a^O.f) be 
successfully*triggered by a shower whose density at the 1 detector 1* p, and 
<x l' y i* “ tho < = oor<,lnate * of the l * h detector. 


We call 

rise) 


R (NO) 


[’ 


sin 20 S(NO) 


m 2W 2 tt 

] ** * / dr / **J * 

0 


d» 


the response (unction of tho array. 

The prograa 1* designed to evaluate thl* (unction for an arbitrary choice of 
the functions a , P , and f, and of detector locations (x , y ) also the Integrations 
can be perforaoi between definite Halts other than 0, 2K , and •» . The Integrations 
are performed by the trapezoidal rule. The density of the net of point* at which the 
Integrand Is evaluated Is adjustable. Experience with hand coaputstlons Indicate* that 
a net containing about 500 points will be required for the desired accuracy In the 
evaluation of P for one set of the parameters and one choice of N and 0. 


The program has been tested with a trial case, and errors are being 
eliminated. 


0. w. Clark 

Nuclear Science Laboratory 


398 A. DIAGONALIZATION OF MATRICES 

In the collective model of tho nucleus, the nucleon* move Independently In 
a potential which represents the major portion of their mutual Interactions. The 
moments of inertia of elllpsoldally deformed nuclei In their ground states are, In 
such a potential, the rigid moments of Inertia corresponding to the rotation of the 
mass distribution about an axis perpendicular to the symmetry axis. The observed 
moments of inertia tend, however, to be less than half this value. The reduction of 
the moments of Inertia has been explained as an effect of that part of the Interaction 
between nucleons not Included In the common potential. 
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-— — .min.•>< mese residual inter 

action, in reducing thm moment, of Inertia, these calculations are being -ad. with a 
greatly simplified model In which the particle, move In an anisotropic harmonic oscllla 
tor potential and tho residual Interactions are represented by an attractive, two-body, 
delta function Interaction. Basis save functions arc formed from single particle har 
monlc oscillator wave functions In Slater determinants and the secular equation Is 
solved approximately (considering only states within the same shell) using the program 
prepared by f. J. Corbato. The elgnefunctIon* are then used to find the moment of 
Inertia of the ground state In the usual way. 

Calculations are being made for two and three particles outside a closed 
shell and for various deformations and strengths of residual interaction. 


It Is expected that work on this problem will be completed early In 1957, 


A. J. Morency 
Physics Department 


399 L,B. DOMAIN WALL MOTION 

In studying the motion of 180° magnetic domain walls In magnetization re 
versal processes In ferrite toroids an equation of tho form 

x + 2d x + f(x) * h 

occurs where x Is the normalized distance the wall moves, 2d Is a damping parameter, 
f(x) Is a non-linear restoring function with maximum and minimum amplitudes of i 1 
snd Is periodic In x with period 2, h is a normalized driving force. The problem is 
to obtain solutions of this equation for various values of d and h. In particular, a 
plot of X or wall velocity vs. time is desired. It Is expected that such plots will 
have the sane general character as the actual voltage response observed In magnetic 
cores during switch processes under the Influence of varying driving currents. 

Most of the time during this quarter was spent In straightening out the 
prograa. One or two good runs have been obtained for certain values of the damping 
constant 8 In the basic equation. Solutions for othor values of 8 and h are to be 
attempted during the coning quarter. 

Programmers working on this problem are J. V. Harrington snd P. Dlneen. 

J. V. Harrington 
Lincoln Laboratory 
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APPENDIX 

1. SYSTEMS ENGINEERING 


WWI RELIABILITY 

The following la the WWI Coaputer Reliability Record for the past quarter: 


Total Computer Operating Time In Hours 1836.5 

Total Time Lost In Hours 35.8 

Percentage Operating Time Usable 98.1 

Mean Tree Time Between failure Incidents (in hours) 23.6 

Total Number of Failure Incidents 84 

Failure Incidents per 24 Hour Day 1,1 

Average Lost Time per Incident In Minutes 25 2 

Average Preventive Maintenance Time per Day in Hours 2,3 


2. VISITORS 


I 


Tours of the Whirlwind I installation Include a showing of the film, 
Making Electrons Count , a computer demonstration, and an Informal discussion 


the major computer 
tailing 237 people 

components. During the past quarter, the folio 1 
visited the computer installation. 

November 

1 

Northeastern University 

27 

November 

6 

Springfield High School 

21 

November 

7 

MIT Courso 2.215, Methods of 
Engineering Analysis 

19 

November 

13 

Naval War College 

13 

November 

15 

Northeastern University 

28 

November 

19 

American Society of Civil Engineers 

113 

December 

17 

MIT Courso 6,25, Machine-Aided 
Analysis'' 

16 


The procedure of holding Open House st the Laboratory on the first of ,h« 
month has continued A total of 56 people attended the 3 Oman Ho... ! “ ! , th 
the quarter, representing -ember, and friend, of th. MIT community Oueen's Collew. 
Simmons College, and the Stromberg-Carlson Company ** College, 
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During the past quarter, there were also 42 Individuals who made brlof 
tours of the computer installation at different times. Represented by those Indi¬ 
viduals were MIT, University of Leyden, Netherlands, U.S. Air Force, U.S. Navy, 
General Electric, Radio Corpcrutlon of America, N.H. Ball Rearing, Inc., Jackson fc 
Moreland, Nippon Electric Company, Blue Cross, Curtiss Wright, American Cyansmld, 
and North American Aviation, 

LET'S TAKE A TRIP 


The weekly Columbia Broadcasting System television program, Lot's Take 
a Trip", was broadcast llvo from the Digital Coaputer Laboratory on Sunday, October 
21 , 1956 at 12 noon, on Channol 12. 

Young Ginger MacManus and Pud Flanagan were led by Sonny Fox on a tour 
through the laboratory. Dr. Francis Bitter, professor of Physics and Associate 
Dean of MIT's School of Science, explained, with the aid of simple experiments, the 
scientific background that led to the development of vacuum tubes, magnetic core 
memory, and other components of the computer. 

Dr. Jay W. Forrester, director of the Laboratory from its inception in 
1951 until July 1956, showed the WHIRLWIND computer in action, explaining the de¬ 
monstration programs which were run by chief operator, Michael Solomlta. Dr. For¬ 
rester also pointed out some of the practical applications of the computer to 
scientific, engineering and Industrial problems. 

Dr. Vannevar Bush, a member of the MIT Corporation, closed the program 
with a few words on the impact that computers have had on science, broadening the 
whole scope of scientific investigation and opening exciting and unlimited pos¬ 
sibilities for the future. 

The show was produced by Stephen Flelschman and directed by Roger Eng¬ 
lander, both of Columbia Broadcasting System. 


3. PUBLICATIONS 


Project Whirlwind technical reports and memoranda are routinely distributed 
to only a restricted group known to have a particular interest in the Project, and to 
ASTIA (Armed Services Technical Information Agency) Document Service Center, Knott 
Building, Dayton, Ohio. Requests for copies of Individual reports should be made to 
ASTIA. 


The following is a list of memoranda published by the Scientific and 
Engineering Computations Oroup during the past quarter. 


No. 


DCL-140 

A Program to Solve the Classical Transportation 
Problem 

J. Dennis 

8-28-56 

DCL-146 

A Proposal for an Algebraic Compllor to bo 

Used on the MIT 704. 

F. Hal Wig 

A. Siegel 

10-16-56 

DCL-148 

List of Available SHARE Material 

F.M. Verxuh 

10-22-56 
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DCL-150 

Simplex Method Subroutine 

K. Bishop 

10-29-56 

DCL-153 

An Argument for Selecting a Certain Size Block 
of Temporary Storage for a Post Mortem Routine 

S. Best 

11-21-56 


Is Explained 



DCL-155 

MIT 704 Post Mortem Routine 

A. Siegel 

11-26-56 


4. ACADEMIC 


1. li u roanct ion 

There were a number ol graduate subject* In automatic computation, numeri¬ 
cal analysis, and data processing offered at MIT during the Fall Term of 1956. Spe¬ 
cifically, the following subjects directly relating to machine computation were of¬ 
fered. 

No. of 


Subject 

Description 

Units 

Year 

Instructor 

Stuc 

2.215 

Methods of Engineering Analysis 

3-9 

G 

Crandall 

17 

6.25 

Machine-Aided Analysis 

3-6 

4 

Arden 

132 

6.567 

Switching Circuits 

3-6 

G 

Caldwell 

52 

15,542 

Management Information Systems 

3-6 

G 

Gregory 

30 

M371 

Operations Research 

3-9 

G 

Wadsworth 

33 

M39 

Methods of Applied Mathematics 

3-9 

G 

Hildebrand 

152 

M411 

Numerical Analysis 

3-2-7 

G 

Hildebrand 

37 


It is apparent from the above list that the subjects are predominantly Graduate A 
subjects—which indicates that the subject is primarily for graduate students. How¬ 
ever, it is perfectly possible for undergraduate students to take these subjects on 
an elective basis. Although the total number of students attending these classes Is 
453, this does not mean that there are that many different people attending these 
subjects since a number of the students are taking several of the subjects. 

2 . 650 Subroutine Developments In the Office of Statistical Services 

Uirlng the Fall Term further developments on Interpretive systems and compi¬ 
lers for the 650 Magnetic Drum Calculator were carried out In the Office of Statisti¬ 
cal Services: 

a) MITSS — MIT Selsctive System 

b) Solution of Linear Systems Using a Modified Elimination Procedure 

c) Evaluation of Hoots of Non-Symmetrlc Matrices Using the 650 Calculator 

d) Polynomial Factorization Research Using the Down-Hill Method. 
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The selective systems have been so srranged to permit a two-pass use of the Symbolic 
Optimal Assemble Program (SOAP) systee on all programs which are written for tbe 650 
Calculator. Details describing the above programs are available In other reports pre 
pared by the Office of Statistical Services. 

3. MIT Computation Center 

Uirlng the Fall Term continued progress has been made In setting up the 
plans for the MIT Computation Center. Specifically, the machine is to be Installed 
In the new Karl T. Compton Laboratory, and the space la rapidly nearing completion. 

It Is expected that the machine will be In operation early in April, 1957. 

During the Fall Term a number of IBM-UIT Research Assistants sere appointed 
as part of this program. They include the following: 


Department 

Name 

Mathematics 

Mr. P. W. Abrahams 

Mechanical Eng. 

Mr. D. M. Baumann 

Civil Eng. 

Mr. C. R. Calladlne 

Physics 

Dr. F. J. Corbat4 

Schl. Ind. Mgment. 

Mr. M. J. Erdel 

Electrical Eng. 

Mr. T. J. Glass 

Chemical Eng. 

Mr. K. Hansen 

Mathematics 

Mr. J. Herahenov 

Physics 

Mr. C. S. Naiaan 

Geology and 

Mr. N. r. Neaa 


Geophysics 
Nuclear Eng. 
Mathematics 


Mr. M. Troost 

Mr. H. J. Welnltschke 


In addition a number of IBM Research Assistants and Associates have been appointed 
at the participating colleges; the following list indicates the present representation: 


Ccllege 

Appointee 

Title 

Boston Coll. 

Dr. W. E. Perrsult 

Mr. J. A. Riley 

Kea. Aaaoc. 
Rea. Aaat. 

Bowdoln 

Dr. R. L. Chlttlm 

Dr. C. E. Huntington 
Dr. E. 0. LaCasce 

Rea. Assoc. 
Kea. Assoc. 

Res. Assoc. 
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PERSONNEL ON THE PROJECT 


College 

Appointee 

Title 

Brown 

Dr. W. Froiborger 

Mr. J. E. Panarelll 

Res. Assoc. 

Res. Asst. 

Conn., U. of 

Mr. M. J. Casarolla 

Rea. Aaat. 

Dartaouth 

Dr. T. 2. Kurtz 

Prof. J. McCarthy 

Res . Assoc. 
lies. Assoc. 

Harvard 

Mr. B. Cohen 

Mr. G. Kraft 

Mr. C. Varsavsky 

Res. Asst. 
Res. Asst. 
Res. Asst. 

New Hamp.. U. of 

Ur. F. J. Loren;cen, Jr. 

Res. Asst. 

Northeastern 

Dr. M. J. Carrabes 
fMr. R. W. Taylor ) 

1Mr. R. F. Trocchlj 

Res. Absoc. 

Item. Asst. 

Tufts 

Dr. K. S. W, Champion 
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A number of programs have been set up for the effective utilization of the 704 
Calculator; these Include a project with the Rockefeller Foundation and a profram 
with the National Science Foundation, 

4. Preparation of Programs for the 704 Computer 

During the Fall Term considerable thought was given to the preparation of 
a suitable post-mortem system for use on the 704 Calculator. This program is cur¬ 
rently being tested out on the General Electric 704 machine at Lynn and it is ex¬ 
pected that it will be available for use iih.nt the 704 arrives in March of 1957, 

A more general algebraic compiler is also receiving attention; however it 
will not be available until later in the academic year. 
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